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Abstract. We examine the characteristic properties of photospheric p- 
modes and sub-photospheric flows of active regions (ARs) observed during 
the period of 26-31 October 2003. Using ring diagram analysis of Dop- 
pler velocity data obtained from the Global Oscillations Network Group 
(GONG), we have found that p-mode parameters evolve with ARs and 
show a strong association with flare activity. Sub-photospheric flows, de- 
rived using inversions of p-modes, show strong twist at the locations of 
ARs, and large variation with flare activity. 
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1. Introduction 



So lar photospheric oscillations havi ng a period of five minutes were first discovered 
byE eighton, Noves & Simonl ([T962). These arise due to p ressure waves (or p-modes) 
which are trapped in d ifferent cavities in the solar interior (luirichl 1 970: Leibac her & Steinl 



I197U; |Peubnerlll975 ). The upper boundary of the cavities lies close to the pho- 



tosphere while the lower one lies at depths depending upon the wavelength of the 
acoustic wave. The global modes of oscillations have lifetimes long enough to travel 
completely around the solar circumfe rence and self interfere without suffering a loss 
of phase coherency greater than tt/2 jHill l 1995). Accurate frequency measurement 



of these modes provides the global properties of solar interior, such as temperature, 
pressure, chemical compositions, etc. But they are insensitive to local characteristics, 
viz., meridional and complex flows of ARs. On the other hand, high degree (I >300) 
p-modes are more sensitive to local properties of ARs. Therefore, photospheric and 
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sub-photospheric proper ties of ARs are studied using loca l oscillations by local heli- 



oseismology techniques ( Lindsey. Braun & Jefferies 19931) . 



Local helioseismology provides a thr ee-dimensi onal view of ARs. Thr ee main 

techn iques are used: ring diagra m analysis (|Hilll 19881). ti me-distance method dDuvall. Jefferies. Harvey & Pomerantzl 
1993) and acoustic holography (Lindsey & Braun 1990). The first one is generaliza- 
tion of the global helioseismology over small area of ARs compared to the whole Sun. 
This is based on the reasonably well understood physics of n ormal modes, while the 
interp retation of other techniques are still in developement (see lGizon. Birch & Spruitl 
2010i and references therein). These methods have been extensively use d earlier 



to stu dy the photospheric p-modes and sub-photospheric flows of ARs (see iMauryal 
2010, and references therein). 



Solar p-modes are believed to be caused by stochastic excitation in the con- 
vection zone but th ey are expec ted to be modified by energetic transients, such as, 



flares and CMEs. IWolfj dl972l) first suggested that a large flare can modify the p- 
mode amplitude by exerting mechanical impulse of the thermal expansion on the 
photosphere. The amplification is expected to be larger in high degree p-modes, 
because they are concentrated near the photosphere, and their typical wavelengths 
approximately matches the scale of the pulse. Subsequently, flare related variations 
in high deg ree p-mode parameters have been observationally reported by several 
researchers d Ambastha. Basu & A ntia 2003; Ambastha. Basu. Antia & Bogartl 12004 



Howe. Komm. Hill. Haber & Hind man 2004; Maurva. Ambastha & TripatfrvOOO' 



ARs are the most important candidates for energetic solar transients. Therefore, 
it is interesting to study their internal structure and dynamics. Local helioseismology 
provides a unique observational tool to determine sub- photospheric flows of ARs by 
inverting the photospheric p-modes (see lMauryall20 1 Oi and references therein). Heli- 
oseismic studies show that sunspots are rather shallow, near-photo spheric phenomena 
( Kosovichev. Duvall & ScherreifeoOOt iBasu. Antia & B ogart 200 ^ and are locations 



of lame scale flows at the photosphere ( Haber et al. 20021: iBraun. Birch & Lindsey 



2004; Zhao & Kosovichev 2004; Komm, Howe, Hill, Gonzalez-Hernandez, Toner & Corbard 
20051: iMaurva & Am bastha 2010a ; b; Maurya, Ambastha & Reddvl bOlhY 



Aim of this study is two fold: First we try to understand the characteristics of 
p-mode parameters at locations of active and quiet regions, and their variations with 
flaring activity. Secondly, we examine sub-photospheric flow variations in ARs from 
non-flaring to flaring phases. We have carried out these studies for the active and quiet 
regions that appeared during 26-31 October 2003. The paper is organized as follow: 
In Section [2) we describe the ARs observed during the aforesaid period. Section [3] 
discusses the observational data and methods of analysis. Section |4] presents results 
of our analysis. Finally, summary and conclusions are provided in Section [5] 
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Figure 1. (a) MDI continuum image of the Sun observed at 11:11:13 UT/28 October 2003. 
NOAA ARs are labeled on their appropriate locations, (b) Integrated GOES X-ray light curve 
in the wavelength range 0.5-4.0A. 



2. Active Regions and Flares 



During the period of 26-31 October 2003, several ARs appeared on the solar disk, but 
the ARs NOAA 10484, 10486, and 10488 dominated due to their complex and large 
size. Amonst these, NOAA 10486 was the largest and most flare productive as it gave 
rise to flares of largest magnitude, viz., X17/4B, X10/2B, of the solar cycle 23. These 
ARs are shown in a MDI continuum image observed at 11:1 1:33 UT/28 October 2003 
(Figure [Ha)). The transient activities for the aforesaid time period are shown by the 
integrated GOES X-ray light curve in the wavelength range 0.5-4. OA in Figure [Hb). 



3. Observational Data and Analysis 



Observational data for the study of ARs for the period 26-31 October 2003 were ob- 
tained from the GONG archive. The p hotospheric p-mode parameters were computed 
using ring diagram analysis (|Hii[l988) of the data cube ( 1 6° x 1 6° x 1 664 m ) . Then, we 
computed horizontal components (u x , u N ) of su b-photospheric flows as a function of 

depth from inversions of p-modes ( Goughl 19851). The vertical component (u z ) of flow 

was c alculated using continuity equation ( Ko mm. Corbard. Durney. Gonzalez Hernandez. Hill. Howe & Toner 
12004 . Maps of kinetic helicity density (KHD), a measure of twist of sub-photospheric 
flows, were then constructed using these three components of flows. 



4. Results and Discussions 



The results obtained from the analysis of ARs for the period 26-31 October 2003 are 
shown in Figure [T]- [3] 
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Figure 2. MDI magnetograms (top) and averaged p-mode (n = - 5) parameters (bottom) 
in the longitude and latitude range ±60° for the period 26-31 October 2003. In the bottom, 
background image corresponds to the frequency shift and vectors represent the horizontal flows. 



4.1 Variations in p-mode Parameters of Active Regions 

Figure Witop) shows the MDI magnetograms in the range ±60° of longitude and lat- 
itude during the period 26-31 October 2003. The bottom panel shows the co-aligned 
maps of the p-mode parameters: frequency shift Sv in the background half-tone im- 
age, and horizontal components of flows Ux and Uy marked by vectors. The ARs 
are labeled in both panels at their appropriate places. The start date/time of the ring- 
diagram data cubes corresponding to the maps are given in the top of the figure. It 
would be interesting to examine the p-mode characteristics as follows: (i) compare the 
flare productive ARs with dormant and quiet regions, and (ii) examine their spatio- 
temporal evolution. For this, we have computed the mode parameters averaged over 
radial orders n = 0-5 in the range ±60° of longitudes and latitudes. 

In Figure ^bottom), vectors represent the horizontal flows and background red 
(blue) colors show the positive (negative) frequency shift in the averaged modes. It 
is evident that there were large, positive frequency shifts at the AR's locations, which 
evolved as the ARs evolved. For instance, NOAA 10488 was initially associated with 
small positive Sv value on 26 October 2003 which increased till 30 October 2003 
as the AR grew in complexity and size. Interestingly, the AR was intensely flare 
productive during this period. 

It is evident that horizontal flows showed a large deviation from the general flow 
pattern around the site of AR NOAA 10486 as seen from the changing direction of 
vectors. Since Ux and Uy are the weighted averages over depth, the flow pattern 
reveals highly sheared flows in the interior of NOAA 10486 as compared to other 
less flare active ARs and the surroundings (see Figure [3] top). These results provide 
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Figure 3. Evolution of (top) magnetic flux and (bottom) kinetic heliciy density during 26-3 1 
October 2003. 



further confirmation on the association of flare productivity of NOAA 10486 with 
p-mode parameters, sub-photospheric flows and their variations. 



4.2 Variations in Sub-Photospheric Flows of Active Regions 

Energetic transients are caused by the changes in magnetic fields which are rooted be- 
neath the photosphere. The sub-photospheric flows result in braiding and intertwining 
of the rising magnetic flux tubes. Therefore, we expect the changes in magnetic field 
configuration to be governed by sub-photospheric flows. In order to investigate this 
issue, we have determined kinetic helicity of sub-photospheric flows corresponding 
to the images shown in Figure [2] 

Kinetic helicity density maps, as a function of depth and latitude for a fixed Car- 
rington longitude 285°, are shown in Figure ^[bottom) during 26-31 October 2003. 
Contour levels are drawn at the 0.5, 2.5, 5, 10, 20, 40, 60, 80% of the absolute max- 
imum of the maps. Top panel shows magnetic flux as a function of latitude corres- 
ponding to the bottom panel. The number of flares in C, M and X class of GOES SXR 
are given in the top panel. The meridian for Carrington longitude 285° passes through 
the center of NOAA 10486 (see Figure 0. Therefore, large magnetic flux (top) and 
kinetic helicity density (bottom) represents the location of NOAA 10486. 
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Figure [3] shows a large twisted flow beneath the ARs as compared to the quiet 
regions. The locations of large KHD (bottom) coincide well with the large magnetic 
flux areas. Magnetic flux in AR NOAA 10486 increased till 29 October 2003 and 
decayed thereafter. This AR was extremely flare productive during 28 and 29 Oc- 
tober 2003 and produced several helioseismic signatures correspondi ng to its major 
flares dDonea & Lindsevll2005l: [KosovichevlEooSlMaurva et al.ll2009 ). On 26 Octo- 
ber 2003, it possessed a large KHD, increasing on 27-28 October 2003, when this AR 
produced the X17/4B flare. After showing a decaying trend, it increased again during 
29 October 2003, when it produced another large X10/2B flare. Thereafter, it simpli- 
fied during 30-31 October 2003. Such systematic variation in internal flows reveals a 
strong role of fluid/magnetic topology contributing in the initiation of large flares. 



5. Summary and Conclusions 



From the study of p-modes and sub-photospheric flows beneath ARs that appeared 
on the solar disk during 26-31 October 2003, we have found the following important 
results: i) Large p-mode frequency shifts occurred at locations of ARs as compared to 
the quiet regions, revealing presence of large surface and sub- surface flows in ARs. ii) 
Sub- surface flows of ARs were comparatively more complex and twisted, iii) Large 
systematic variations in internal flow parameters were found associated with flaring 
activity. 
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